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The statistical mechanics of a two-state Ising spin-glass model with finite random connectivity, 
in which each site is connected to a finite number of other sites, is extended in this work within the 
replica technique to study the phase transitions in the three-state Ghatak-Sherrington (or random 
Blume-Capel) model of a spin glass with a crystal field term. The replica symmetry ansatz for 
the order function is expressed in terms of a two-dimensional effective-field distribution which is 
determined numerically by means of a population dynamics procedure. Phase diagrams are obtained 
exhibiting phase boundaries which have a reentrance with both a continuous and a genuine first- 
order transition with a discontinuity in the entropy. This may be seen as "inverse freezing", which 
has been studied extensively lately, as a process either with or without exchange of latent heat. 
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INTRODUCTION 



Numerous problems of frustrated, disordered systems, with extensive connectivity in which each site is linked to a 
macroscopic number of other sites (l| have been studied in the past by means of the replica technique in mean-field 
theory @, S|- The technique has been extended to systems with finite random connectivity and binary units (spins) 
in states a = ±1, in which each site is linked to a finite number of other sites, in areas like error correcting codes 
[US; HI; spin glasses neural networks [l2-14| and small- world lattices [la]. The fact that mean- field theory is 



exact only for infinite-range interactions or infinite-dimensional systems has been a challenge for the understanding of 
the behavior of more realistic disordered systems, and a study of the effects of finite connectivity even in mean-field 
theory could be a useful improvement. 

Three-state spin models of states a = 0,±1 with random bonds and finite connectivity could be of interest to 
condensed matter physics in view of the phase transitions that already appear, within mean- field theory, in the 
random Blume-Emery-Griffith-Capel (BEGC) model with full connectivity between the spins [l6, 17 1 and it should 



also be of interest for information processing in neural networks [181 ] . The simplest case of a fully connected random 
BEGC model is the three-state sp in-g lass model of Ghatak and Sherrington (GS) flil. [20|. which is a Blume-Capel 
(BC) model with random bonds [2l|. The GS model with infinite range interactions exhibits both a continuous 
transition at high T and a genuine thermodynamic first-order transition below a tricritical point between a spin-glass 
and a paramagnetic phase. The first-order transition appears in a reentrant part of the phase boundary and it may 
describe "inverse freezing" with an exchange of latent heat. This is a reversible transition, say between a paramagnetic 
(P) phase and a spin-glass (SG) phase in which the entropy of the P phase below the transition is smaller than the 
entropy of the SG phase and there has been a recent revival of interest in these transitions (22 , 23 1 which could explain 



the behavior of colloidal and polymeric systems, among others (see ref. [23j for a recent summary of realizations). 

A result of the GS model, either in mean-field theory with full connectivity or obtained by means of numerical 
simulations for nearest-neighbor interactions in three dimensions, is that the tricritical point separating the continuous 
and the first-order transition is either above or at the point of reentrance on the phase boundary implying that inverse 
freezing appears only as a first-order transition flil . Il7l . |23| . 

General results in the form of complete phase diagrams of stationary states for the GS (or random BC) model with 
finite connectivity are still missing and the purpose of this paper is to make progress in that direction by means of an 
analytic study extending the SG replica technique, with the replica-symmetry (RS) ansatz, developed for disordered 
Ising systems with finite random connectivity and binary spins [T^ . Typical questions one would like to answer is 
how the nature of the transition changes with the connectivity, if there is reentrance behavior of the phase boundaries 
and if this takes place either as a continuous or as a first-order transition at finite temperature. 

The procedure on which the present work is based introduces an order- function 0] , in place of an infinite number of 
order parameters Q , and uses a representation in terms of a weighted probability of alignment of the spins involving 
a distribution of a two-component local field. One of the components is associated to a linear form in the spins, 
while the other component is associated to the crystal field term. A population dynamics technique is used to solve 
numerically a self-consistency equation for the distribution of the local field that yields the relevant physical order 
parameters. Phase diagrams are then obtained which exhibit either a continuous transition between a P and a SG 
phase, for low connectivity, or both a continuous and a first-order phase transition that merge at a tricritical point 
with reentrance behavior characteristic of inverse freezing between those phases, for increasing finite connectivity. 

The outline of the paper is the following. In section 2 we present the GS model and carry out the replica procedure 
that yields the statistical mechanics for the three-state system. In section 3 we present the results for the distribution 
of local fields, the relevant order parameters and the free energy and discuss the phase diagrams. We conclude in 
section 4 with a summary and remarks. 



THE MODEL AND REPLICA PROCEDURE 

We consider a system of N interacting three-state Ising spins ui E { — 1, 0, 1}, i = 1 . . . N, described by the Hamil- 
tonian 



H = --^c ij J l] a l a j + D^a 2 . (1) 
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The random variable c^- € {0, 1} indicates whether there is a connection (c^ = 1) or not (cy = 0) between a pair of 
spins and it takes different values for different pairs of spins, according to the distribution 



(2) 



where c (the connectivity) is the average number of connections per spin which is assumed to remain finite in the 
thermodynamic limit N — > oo, such that c/N — > 0. Thus, the sites are connected according to a Poisson distribution, 
and one makes use of this limit in deriving the statistical mechanics of the system. There is a set of infinite-range 
interactions {Jij} that will be assumed to be independent, identically distributed, random variables drawn from a 
distribution p(Jij), to be specified below, and averages over that distribution of a quantity g(Jij) will be denoted by 
(g(J)) j- The quadratic form in the spins favors the population of the zero state, if D > 0, or the states ±1, if D < 0. 
If D is sufficiently large and negative one retrieves the binary Ising spin-glass model with finite connectivity and spins 
oi 6 { — 1,1}, which is a particular case of the recently studied small world spin glasses when the nearest neighbor 
interaction along the ring is set to zero (l5| . 

Assuming thermal equilibrium at an inverse temperature (3 = 1/T, the disorder-averaged free energy per spin is 
calculated in the replica procedure as 



f(fi) = - lim lim - log(Z n ) 



where 



E 



-8H 



(4) 



is the partition function and the brackets stand for the disorder average. In the small c/N limit, the disorder-average 
replicated partition function becomes, to leading order in N, 



(Z n ) 



E 
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7 > - 1 



(5) 



where a = 1, . . . , n denotes the replica index. Since the connectivity c is finite, one cannot expand the inner exponential 
and follow the standard infinite-connectivity calculation. Instead, to extract the variables under summation from the 
inner exponential, one introduces the identity 



CT a=l 



where cr and cr^ are n-component vectors representing replica states and Scrcr, 
Thus, we write 



(6) 



1 if cr — cr, and zero otherwise. 



(Z»)= ]T exp 
a 1 a n 



i& (TT 



(7) 



and introduce the order function P(cr) = (1/N) J^- £<T<T i5 which represents the fraction of sites with the replica 
configuration cr, through the identity 



1 = f]ldP(*)dP(cT)eZ° p W( p W-^ 5 ™*) 



(8) 



where P is an auxiliary density. Performing the trace and changing P to NP the integral 



Z n ) = I Y[dP(cr)dP((T)cxpN 
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can be evaluated by the saddle-point method in the large- TV limit, leading to the extremum 



f(B) = - lim 4-Extr 

K ' n->0 fin 



Y,P(*)P(<T)+logY,e-^ a *l-P«T) 
(J (J 

+ ^P W P(r)(c^^--l 



over the densities {P(er), P(er)} for the free energy per site. The saddle-point equations become 



P(<r) 



e-P D T, a °l-Hi 



and 



(10) 



(11) 
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Eliminating P(er) in Eq. (jTTJ) by means of this expression we obtain the self-consistency relationship 

exp (-02? £ a ^ + c£ T (e^ J S«, _ x ) ) 



P(<r) = 



E CT exp (-/3£> £ Q + c£ T P( T ) (e^ J S„ *«t„ _ ^ 



and the free energy becomes 



lim — 



|^P( CT )P(r)(e^^— -1 
err 

log 2 exp L^^^ + c^P(r)(e^^— -l) ) 



(13) 



(14) 



with the order function given by the saddle-point equation, indicated by the asterisk. 

Our search for solutions of Eq. (fT3")) will be restricted to the replica symmetry (RS) ansatz, This means 

that P(er) should remain invariant under replica permutations and, consequently, for three-state spins it should only 
depend on the summations £ Q a a and £ Q (er a ) , with their correspondin g w eights. Thus, in extension of the RS 
ansatz for the order function in finite-connectivity two-state Ising models [la |. we assume that for the three-state 
model 



P(cr) 



dh dbW(h,b) 



[2e-P b coMfih) + 1] 



(15) 



for any real n, where h and b are the two components of the local field and W(h, b) is a density which has to be 
determined self-consistently. Since the normalization factors, both here and in Eq. (|13p . become one in the limit n — > 
we may leave them aside. Expanding the second exponential in Eq. lflU)) and using the RS ansatz we have 



P(<r) 



k\ 



E 

k=0 



Y[dh t dbiW (hi,h)dJip(Ji) 

^log^e /3/tiT °- /3&! ( T °) 2+ -' 3J,a ° T ° 



k 

JJexp 



(16) 



where p (J) is the probability distribution for the interaction, assumed to be continuous. Summing over r l a and using 
the identity J2 a ^<r<r a = 1, in order to extract the appropriate dependence on £ Q <j a and £ Q (a a ) , we obtain 

e c 



— — C K f 

P( <7 ) = e-^E-^^l-£- / WdhdhWihub^dJwiJi 

k=0 ' 1 = 1 
k 



I cx p J2 5 °°° lo § I 2 cosh + 0-Wc) e ~^ bl + 1] 



(17) 
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For three-state spins we use the representation 

S aaa = l-(J 2 -cjl + aa a /2 + 'icj 2 al/2 (18) 

to do the sum over a and with the RS ansatz in the left hand side of the equation we obtain the self-consistency 
relationship for W(h,b), 

oo _ c - k 

W(h,b) = ^^- HdhidbtWihiAWipiJi) (19) 

k=0 ' J 1 = 1 

8(h-~J2t(htA,Jl))s(b-D + ~J2i>M,Ji)) , (20) 



where 



and 



H 1=1 / \ ' 1=1 



1 2ai -I- e /3b 



, /, . n 1 , (2 fll +c^)(2a_ 1+ e^) 

^ (h, b, J) = - log —2 , (22) 



2 



[2a +e^] 



in which a a = cosh(/3/i + a(3J/c) . 

To determine the density W (h, b) we proceed numerically by means of population dynamics of a large number of 
fields updated as follows [ll|, for each value of (3 and D. First a number k is chosen from a Poisson distribution of 
mean c. Then, cells (hi,bi) and couplings Ji with I running from 1 to are selected at random from the population 
and the summations in the delta functions are calculated. Next, one selects at random a new cell (h, b) from the 
population and sets 



k 

= ^^(/>(hi,b h Jl) , (23) 



i=i 



r i=i 

continuing the procedure until it converges to a limiting W{h, b). 
Knowledge of W(h, b) allows to determine the magnetization 



1 k 

b = D--^(h h b l ,J l ) (24) 



m = / dhdbW(h,b) (a) , (25) 



the spin-glass order parameter 



and the additional parameter 



q= I dhdbW(h,b)(a) 2 (26) 



dhdbW(h,b)(a 2 ) , (27) 



for the three-state model, where (a) and {V 2 } depend on h and b as 

sinh(/3/i) 



3h(/3/i) + e/ 3b /2 



(28) 



(i 




(a) (b) 

FIG. 1. Local field distribution W(h,b) for c = 6.0, D = 0.35 and T = 0.1, with initial conditions (A) and (B) specified in the 
text in (a) and (b), respectively. 



and 



<- 2 > 



cosh(/3/i) 



cosh(/3/i) + e^/2 



(29) 



As usual, m 7^ indicates magnetic ordering, while m = and q > corresponds to spin-glass ordering. The 
additional parameter r is zero only when all spins are in the a = local state. 
The free energy has energetic and entropic contributions given by 

/ 1 ( ( g) = .± f dhdti dbdb' dJ P(J)W(h,b)W(h',b') 
J 

x /log 1/4 + c-^ (6+6 ' ) (^cosh(/3/i)cosh(^')cosh(/3 J/c) 
+ siah(0h) ainh(j3h') sinh(^ J/c)) +c~ pb cosh(/3/i) /2 + e"' 36 ' cosh(/3/i') /2 



and 



-log 



e"' 36 cos 



hOSfc) + 1/2) (c- /3b ' cosh(/3/i') + I/ 2 )] } 



1 c fc e 



fc=0 
X lOE 



A-! 



Jdhi db e dJ f W(h e ,b e )P(J e ) 



I + 2exp(^2i> (he,b e , Jj,) - pD \ cosh( ^ (h,, 6/, J t ) 



(30) 



(31) 



respectively, within the RS ansatz. 



RESULTS 



In what follows we assume that the couplings Jij are independent Gaussian random variables with zero mean and 
unit variance. In the binary Ising SG model, there is a one component local field h and a density W(h) = 5(h) in 
the P phase. The critical temperature for a bifurcation to a SG solution can be found in that case by means of an 
expansion of the full W(h) around 6(h). (See ref. [IBJ for details.) In the case of the present three-state SG model 
with finite connectivity, there is no simple form for the density of the two-component local field W(h,b), even in the 
P phase, and in order to obtain the thermodynamic properties one has to use a population dynamics procedure to 
calculate explicitly that density. Since the mean of Jij is zero, there is no long-range order and we expect that W(h,b) 
is symmetrically distributed around the h — axis. Depending on the values of r and q, the system can be found 
either in a SG or in a P phase. 
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D 



FIG. 2. Order parameters r and q and free energy per site / as functions of D, for c = 8.0, exhibiting a continuous transition 
at T = 0.07 (light solid lines) and a first-order transition at T = 0.04 (heavy solid lines for SG states and dashed lines for P 
states) discussed in the text. The arrows indicate points on the spinodals. 

The implementation of population dynamics requires an initial guess for W(h,b), to be constructed as follows. A 
finite surface on the space (h,b) is divided into n x n cells and this set of n 2 cells is populated with N vector fields. 
In this work, the initial population was distributed in two distinct ways: (A) the two components of each field were 
randomly distributed between -0.5 and 0.5 and (B) all the fields had h = and b randomly distributed between 
D — 1 and D. For each set of parameters c, D and T the population dynamics is allowed to run until a stationary 
distribution W(h,b) is reached. 

Illustrative examples of stationary distributions for c = 6, D = 0.35 and T = 0.1, computed with n — 128 and 
N = 40 000 are shown in Fig. [TJ Starting from the initial condition (A), the population dynamics converges to 
the two-dimensional distribution, symmetric around h = 0, shown in Fig. [T](a). This field distribution leads to a 
spin-glass solution with r > and q > 0, according to Eqs. (|2T| and ([2"fJ|) . Starting with the initial condition (B), 
W(0,b) converges to the one-dimensional distribution shown in Fig. [1] (b). Here, only fields along the h = axis 
remain populated, and there is a long tail along the b axis. This field distribution leads to a paramagnetic phase with 
r > and q = 0. Thus, for the chosen set of parameters, the initial conditions (A) and (B) are within the basins of 
attraction of the SG and P solutions, respectively. 

Specific examples of results for the dependence of the order parameters with D near the transition between the SG 
and the P phase obtained from Eqs. (|26[) and (|27[) . that were used to construct the phase boundaries are shown in 
Fig. H]for connectivity c = 8.0 and two typical situations, one of the higher and the other one of the lower temperature 
behavior in which T = 0.07 and T = 0.04, respectively. The corresponding free energy for both situations, obtained 
by means of Eqs. (30) and (|3"Tj) . is shown in the lower panel of that figure. Each point on the curves is a result of an 
average over ten runs of the population dynamics procedure, in order to smooth out fluctuations. In the presence of 
multiple solutions for the parameters, we follow previous works on the fully connected model and choose the largest 
one in magnitude, that yields the lowest free energy which can be smoothly continued from one phase to the other 
with a change in D. The curves in Fig. [5] are rather close to those obtained already with a lower resolution, say 
n = 64 and N = 10 000, indicating that the dynamics is near convergence. To test this point we found that further 
results for r, q and / obtained with n = 256 and TV = 160 000 at T = 0.04, which is the case where the transition is 
first-order, are almost indistinguishable from those in Fig. [2 

As D increases for T = 0.07, q and r decrease continuously and the transition from the SG to the P phase takes 
place at q — and small r. When T = 0.04, instead, there is a genuine first-order transition signalled by a continuous 
free energy with a discontinuity of the entropy and the appearance of a pair of spinodals, one for states (in solid 
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FIG. 3. (a) The T vs. D phase diagram for connectivity c = 2.0, 4.0, 5.6, 6.0 and 8.0 from right to left. Solid (dotted) lines 
denote continuous (discontinuous) transitions between a spin glass (SG) phase at left and a paramagnetic (P) phase at right of 
each curve, (b) Detail of the phase boundary for c = 8 including the spinodals of the P and SG phases on the left and right of 
the transition in light solid lines, respectively. 

lines) attained from the spin-glass phase and the other one for states (in dashed lines) reached from the paramagnetic 
phase with increasing or decreasing values of D, respectively. The arrows indicate the points on the spinodals for 
that value of T, and the transition between the phases appears at the crossing of the free energies, where D = 0.328. 
The first-order transition persists for somewhat higher T with the spinodals becoming closer up to the merging with 
the continuous transition at a tricritical point given by T ~ 0.059 and D ~ 0.33 for c = 8.0, and similar results are 
obtained for other values of c > c* where c* is between 5.6 and 6.0, whereas only a continuous transition appears for 
smaller c at all T. It is interesting to note that the changeover from a continuous to a discontinuous transition at low 
T is preceded by a discontinuous transition between two SG states within the SG phase for c* ~ 5.6. For c slightly 
larger this transition merges with the continuous transition that separates the two phases. 

A global picture of the transitions for finite connectivity is given by the T vs. D phase diagram for several values of 
c < 8.0 shown in Fig. [3J The transition is continuous for low c and all T with the order parameters going continuously 
to zero with increasing values of D, and there is a reentrance for c > 3.5. On the other hand, a first-order transition 
appears with increasing c between 5.6 and 6.0 at low but finite T that starts at a tricritical point on the reentrance 
of the continuous phase boundary. The reentrance separates the P phase at low T from the SG phase at higher T. 
Note that the tricritical point appears at larger values of T with increasing c, and this has been checked by further 
calculations for c = 16. We also found that, for all the values of c larger than ~ 2.5 considered in this work, the phase 
boundaries move to the left with increasing c towards lower values of D. Indeed, there is a monotonic decrease of 
Tq(c), which is the critical T at D = 0, with increasing c and we found numerically that To ~ J/y/c for large c, which 
is the same behavior as that in the two-state spin-glass model with finite connectivity where the transition is always 
continuous (l5| . 

The reentrance of the phase boundaries is a manifestation of inverse freezing which appears either as a genuine 
first-order transition with a discontinuity of the entropy or as a continuous transition with only a gradual change in 
the entropy from one phase to the other. In order to show that the discontinuity of the entropy across the first-order 
transition in the GS model with finite connectivity goes in the right direction to account for inverse freezing we show 
in Fig. U] (a) the temperature dependence of the free energy for c = 8 and D = 0.325. The free energy is continuous 
and follows the lower curve when there are two solutions and, indeed, the entropy of the P phase below the first-order 
transition (which is slightly positive for that value of D) is smaller than the entropy of the higher temperature SG 
phase, with a discontinuity of the entropy at the transition. This means that in heating the system the disordered 
paramagnet "freezes" into the amorphous SG with an exchange of latent heat. For comparison, we show in Fig. 0] 
(b) the free energy, for D = 0.35, with no discontinuity in the entropy across the continuous transition. 

The phase diagram for finite c obtained in the present work within the RS scheme differs from the T vs. D phase 
diagram for the fully connected GS model in that, apparently, there is no reentrance in the latter within that scheme 
but there is a reentrance within full replica symmetry breaking. The first-order transition line appears in that case 
all along the reentrance of the phase boundary down to T = and it is even slightly continued above into the normal 
phase boundary. Thus, there is no continuous transition on the reentrance in that case. On the other hand there is, 
apparently, a further reentrance on the first-order transition at low T [13]. 
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FIG. 4. Temperature dependence of the free energy for c = 8 across the first-order transition at D = 0.325 in (a) and the 
continuous transition at D — 0.35 in (b). The dashed left (right) line indicates a point on the SG (P) spinodals, respectively, 
and the vertical solid lines indicate the transition with an increase in entropy from the P to the SG phase in (a) and with no 
increase in (b). The normal transition between the SG and the P phase, not shown in the figure, appears at higher T. 



The kind of behavior we find here for finite c is similar to the phase diagram of a somewhat different model of a 
three-state Ising spin glass with full connectivity in which the degeneracy of the active spins [a = 1 or a = — 1) is 
larger than the degeneracy of the inactive state (<r = 0), in contrast to the GS model where these degeneracies are 
the same. Indeed, the T vs. D phase boundary in that case also displays both a first-order transition at low T and a 
continuous transition at higher T on the reentrance of the phase boundary, even within the RS scheme [22| . 



SUMMARY AND CONCLUSIONS 



The statistical mechanics of the binary Ising SG model with finite random connectivity, within mean-field theory, 
has been extended in this work to study the three-state Ising SG model with crystal-field effects of Ghatak and 
Sherrington, which is a random-bond Blume-Capel model. Renewed interest in this model is due to the recent 
discovery of numerous physical realizations within condensed-matter physics that exhibit inverse transitions (see ref. 
[23| for a recent review), and finite connectivity could be of use in order to study the behavior of such systems with 
effective interactions between infinite and short range or between infinite dimensions and three-dimensions. 

The replica method for disordered spin systems with finite connectivity has been extended here within the replica 
symmetry (RS) ansatz. The RS ansatz for the order function that describes the behavior in the present model 
with finite connectivity introduces a distribution of a two-component effective field which plays a crucial role in 
determining the relevant order parameters of the model. An iterative self-consistency relation is derived for the local 
field distribution, an analytic solution of which does not seem possible in the presence of a first-order transition, 
and we use instead a population dynamics numerical procedure. Phase diagrams were then obtained in order to 
investigate the effects of finite connectivity and anisotropy parameter on the transitions between a SG and a P phase, 
with emphasis on reentrance behavior. The latter is a necessary feature in order to describe inverse transitions, in 
particular inverse freezing, in view of the two phases favored by the present model. 

The main result of this paper is the T vs. D phase diagram, for various values of the connectivity, which exhibits the 
different phase transitions that may appear at low temperature. For low connectivity the transition is a continuous one 
down to T = 0, with a reentrance at low T for higher connectivity exhibiting both a continuous and a discontinuous 
transition. We showed that the discontinuous transition is a first-order transition with a discontinuity of the entropy 
at the transition in which the entropy of the SG phase is larger than the entropy of the P phase below the transition, 
which is a sign of inverse freezing in the system. 

The results obtained so far within the RS ansatz could be changed by replica-symmetry breaking (RSB) and an 
extension in this direction could be interesting although not so easy to carry out in the present case due to the time 
consuming evaluation of the distribution of the two-dimensional effective field. However, if the trend of the results for 
the GS model with full connectivity are taken as a guide, one would expect an enhancement of the reentrance on the 
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first-order transition for finite connectivity within RSB, without much change of the continuous transition. 

It may also be interesting to find out to which extent the specific form of the distribution of random bonds makes 
a difference. In place of a Gaussian with mean zero and unit variance that we used in this work, one could have 
a Gaussian with mean Jo > 0. We found that the phase diagram is somewhat changed in that case. Instead, one 
could consider a general bimodal distribution between ferromagnetic and antiferromagnetic interactions which allows 
to predict when the effects of frustration may become more relevant before going into a RSB calculation, an argument 
that has been used before [la ]. This, and related issues, are currently being studied. 

The procedure extended in this work can be applied to study the effects of finite connectivity in other multi- 
state spin models like the random BEGC model or the degenerate spin model of Schupper and Shnerb either with 
random or uniform interactions (l7l . |22| . It can also be applied to attractor neural networks and the study of the 
retrieval performance in a three-state network with finite connectivity and a Hebbian learning rule is currently being 
investigated [24j |. 
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